set more off

use dataset, clear

local Z pcteligible2013

keep if abs(x) < 100

su `Z', det
gen Z = `Z' >=r(p50)
replace Z = . if mi(`Z')

label define zl 1 "High eligibility" 0 "Low eligibility", replace
label values Z zl

#delimit

gr tw 
	(sc dpctui20142013 x, msize(tiny) col(gs13))
	(lpolyci dpctui20142013 x if x < 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
	(lpolyci dpctui20142013 x if x > 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
	,
		ylab(-40(10)10, angle(horiz) nogrid)
		plotregion(style(none))
		subtitle(, fcol(none) lwid(none))
		xline(0)
		ylab(, angle(horiz))
		xtitle("Distance from nearest state border")
		legend(off)
		ytitle("Change in pct. uninsured")
		title("2014-2013")
		legend(off)
		nodraw
		name(g1, replace)
		;		
	
#delimit cr

#delimit

gr tw 
	(sc dpctui20152013 x, msize(tiny) col(gs13))
	(lpolyci dpctui20152013 x if x < 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
	(lpolyci dpctui20152013 x if x > 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
	,
		ylab(-40(10)10, angle(horiz) nogrid)
		plotregion(style(none))
		subtitle(, fcol(none) lwid(none))
		xline(0)
		ylab(, angle(horiz))
		xtitle("Distance from nearest state border")
		legend(off)
		ytitle("Change in pct. uninsured")
		title("2015-2013")
		legend(off)
		nodraw
		name(g2, replace)
		;
	
#delimit cr

gr combine g1 g2, rows(1) title("Change in uninsured") name(G1, replace) 

gr export "_output/figureA7.pdf", replace
gr close
